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We report results of our ongoing nonperturbative numerical study of a classical ef- 
fective theory describing low-x partons in the central region of a heavy-ion collision. 
In particular, we give estimates of the initial transverse energies and multiplicities 
for a wide range of collision regimes, including those at RHIC and at LHC. 

In the central-rapidity region of heavy-ion coUisions at RHIC and at LHC 
a combination of very high center-of-mass energy with a very large number of 
participating valence quarks is expected to lead to a novel regime of QCD, one 
characterized by a very high parton density. This regime may not be amenable 
to analysis by conventional methods, such as multiple scattering fsi classical 
cascade descriptions, which ignore the coherence of the gluon fieldEl. 

This coherence emerges naturally from the_classical effective field theory 
approach of McLerran and Venugopalan (MV) □. If the parton density in the 
colliding nuclei is high at small x, classical methods are valid. It has been shown 
recently that a RG-improved generalization of this effective action reproduces 
several key results in small-x QCDu. 

The model can be summarized as follows. The high-x and the low-x par- 
tons are treated separately. The former corresponds to valence quarks and 
hard sea partons. These high-x partons are considered recoilless sources of 
color charge. For a large Lorentz-contracted nucleus, this results in a static 
Gaussian distribution of their color charge density p in the transverse plane: 



The variance of the color charge distribution is the only dimensional pa- 
rameter of the model, apart from the linear size L of the nucleus. For central 
impact parameters, can be estimated in terms of single- nucleon structure 
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functionsu. It is assumed, in addition, that the nucleus is infinitely thin in the 
longitudinal direction. Under this simplifying assumption the resulting gauge 
fields are boost-invariant. 

The small x fields are then described by the classical Yang-Mills equations 

D^,F^, = J, (1) 

with the random sources on the two light cones: — ^ ^v,±S{x^)pi,2{i't)- 
The two signs correspond to two possible directions of motior^along the beam 
axis z. As shown by Kovner, McLerran and Weigert (KMW)El, low- a; fields in 
the central region of the collision obey sourceless Yang-Mills equations (this 
region is in the forward light cone of both nuclei) with the initial conditions in 
the At =0 gauge given by 



A' = + A. 
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A^ = 

2 



(2) 



Here the pure gauge fields A\ ^ are solutions of (|^) for each of the two nuclei 
in the absence of the other nucleus. 
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Figure 1: Transverse-plane energy density per unit rapidity versus proper time for the values 
5.66 (diamonds), 35.36 (plusses), and 297 (squares) of As-L.Both the energy density and the 
proper time are expressed in units of As . The solid lines are fits of the data to the form 
a -\- /3exp(— 7t). 



In order to obtain the resulting giuon field configuration at late proper 
times, one needs to solve (|l|) with the initial condition (||). Since the latter 
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depends on the random color source, averages over realizations of the source 
must be performed. KMW showed that in perturbation theory the gluon 
number distribution by transverse momentum (per unit rapidity) suffers from 
an infrared divergence and argued that the distribution must have the form 



for k± 3> As. The log term clearly indicates that the perturbative description 
breaks down for k± ^ A^. 

A reliable way to go beyond perturbation theory is to re-formulate the 
EFT on a lattice by discretizing the transverse plane. The resulting lattice 
theory can then be solved numerically. We shall not dwell herejpn the details 
of the lattice formulation, which is described in detail in Ref. eI. Keeping in 
mind that A^ and the linear size L of the nucleus are the only physically in- 
teresting dimensional parameters of the modelLI, we can write any dimensional 
quantity q as Af/g(AsL), where d is the dimension of q. All the non-trivial 
physical information is contained in the dimensionless function fq{KsL). We 
can estimate the values of the product AgL which correspond to key collider 
experiments. Assuming Au-Au collisions, we take L = 11.6 fm (for a square 
nucleus!) and estimate the standard deviation A^ to be 2 GeV for RHIC and 
4 GeV for LHCa. Also, we have approximately g ^ 2 for energies of interest. 
The rough estimate is then A^L « 120 for RHIC and A^i « 240 for LHC. 

We now proceed to describe the results of our numerical study so far. Our 
simulations were performed for the SU(2) gauge group, in order to keep the 
computational costs low. 



Table 1: The function /b, i-e., the energy per unit transverse area per unit rapidity, expressed 
in units of A^, tabulated vs AgL. 

We first compute the energy per unit transverse area per unit rapidity, 
deposited in the central region by the colliding nuclei. As Figure ^ illustrates, 
this quantity tends to a constant at late proper times. We find that at r — > oo 
the energy density in units of As depends on the dimensionless parameter A^L 
as described in Table 1. Note the very slow variation of this dimensionless 
function in the entire range of AgL values, which includes both our RHIC and 
LHC estimates. Using this plot, and assuming, in accordance with Ref. the 
{N^ — l)/Nc dependence of the energy on the number of colors Nc, we arrive 




(3) 



AsL 

fE 



17.68 35.36 70.7 106.06 148.49 212.13 296.98 
.323(4) .208(4) .200(5) .211(1) .232(1) .234(2) .257(5) 
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at the values of 2700 GeV and of 25000 GeV_for the transverse energy per unit 
rapidity at RHIC and at LHC, respectivelyB. 
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Figure 2: n(k^) = dN / L'^ /d?k± as a function of the gluon momentum k±_ for KsL = 297 
(diamonds). The sohd hne is a fit to Bose-Einstein distribution, corresponding to a gas 
of free particles at the inverse temperature f) = (1.28 it 0.01)As and with the mass m = 
(0.092 ±0.003)As. 

The number and distribution of produced gluons are of considerable in- 
terest as initial data for possible evolution of the gluon gas towards thermal 
equilibriumlj'tjEil. Strictly speaking, a particle number is only well-defined in 
a free-field theory, and there is no unique extension of this notion to a gen- 
eral interacting case. For this reason we use two different generalizations of 
the particle number to an interacting theory, each having the correct free-field 
limit. We verify that the two definitions agree in the |5iieak-coupling regime 
corresponding to late proper times in the central regionllj. 

Our first definition is straightforward. We impose the Coulomb gauge con- 
dition in the transverse plane: V_l • A±^ = and and determine the momentum 
components of the resulting field configuration. Our second definition is based 
on the behavior of a free-field theory under relaxation. Consider a simple 
relaxation equation for a field in real space, 

dtcp{x) = -dV/d^{x), (4) 

where t is the relaxation time (not to be confused with real or proper time) and 
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V is the potential. It is then easy to derive the following integral expression 
for the total particle number of a free-field system: 




where V{t) is the potential for the relaxed field. Now (P can be solved numer- 
ically for interacting fields. Subsequently, V{t) can be determined, and N can 
be computed by numerical integration. This technique presently only permits 
determination of the total particle number. 

We now present our results using both the techniques discussed. Our 
results for the number distribution, computed in Coulomb gauge, are as follows. 
We have verified that for large k± our numerically obtained multiplicity agrees 
with the lattice analogue of the perturbative expression . At smaller k± , the 
distribution softens and converges to a constant value, unlike its perturbative 
counterpart. Notably, this qualitative change of the distribution occurs at 
A:_L ~ As. We tried to quantify this non-perturbative behavior by fitting the 
distribution to a variety of physically motivated functional forms. Surprisingly, 
we find that the shape of the distribution is closely reproduced by the Bose- 
Einstein form n(k) = A/(exp(/3wfc) — 1), with the inverse temperature (3 of 
the order of 1 in units of As, and with cok corresponding to a free massive 
dispersion relation, with the mass of the order of 0.1 As. This is an unexpected 
result for a purely classical theory, whose meaning, beyond providing us with 
a convenient parametrization, is not yet clear. 

Table 2 summarizes our findings for the total gluon number, written as 
N = (As£/g)^/Ar(AsL). Note the good agreement between the two methods 
used to determine N (a more detailed analysis shows that the small difference 
between the two is likely an artifact of the Coulomb gauge). Note that, simi- 
larly to fsi-^sL), /jv(AsL) is a slowly-varying function. Based on Table 2, and 
making a naive extrapolation from the SU(2) to the SU(3) gauge group, we 
estimate the initial gluon number per unit rapidity as ^ 950 for central Au-Au 
collisions at RHIC energies, and as ~ 4300 for central Au-Au collisions at LHC. 
Let us again emphasize that, even though our figure for RHIC is not far from 
the recent experimental result, both the distribution and the total multiplicity 
are likely to be changed by further evolution of the field configuration. 
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Table 2: Values of /jv vs A^L for the two definitions of the total particle number. 
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At the current stage of our project we are only able to make qualitative 
predictions. One obvious way to improve the accuracy is by replacing the SU(2) 
color group by the physical SU(3) one. Another is by allowing deviations from 
the strict boost invariance. These issues will be addressed in the future. 
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